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Abstract. In string theory the couphng "constants" appearing in the low-energy 
. . . effective Lagrangian are determined by the vacuum expectation values of some (a 

' priori) massless scalar fields (dilaton, moduli). This naturally leads one to expect 

, a correlated variation of all the coupling constants, and an associated violation of 

' the equivalence principle. We review some string-inspired theoretical models which 

incorporate such a spacetime variation of coupling constants while remaining natu- 
rally compatible both with phenomenological constraints coming from geochemical 
data (Oklo; Rhenium decay) and with present equivalence principle tests. Barring a 
very unnatural fine-tuning of parameters, a variation of the fine-structure constant 
QQ ■ as large as that recently "observed" by Webb et al. in quasar absorption spectra 

appears to be incompatible with these phenomenological constraints. Independently 
of any model, it is emphasized that the best experimental probe of varying constants 
are high-precision tests of the universality of free fall, such as MICROSCOPE and 
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1. Introduction 



String theory is the only known framework which reconciles pertur- 
bative gravity with quantum mechanics. Though string theory is not 
fully defined (beyond the perturbative level), and though there is no 
^! clear understanding of the way the Standard Model of particle physics 

eventually fits within this theory, it seems interesting to take seriously 
some of the hints suggested by the presently understood parts of string 
^ , theory. One of these hints is that the coupling "constants" appearing 

■ in the low-energy Lagrangian are determined by the vacuum expec- 

tation values (VEV) of some a priori massless scalar fields: dilaton 
and moduli. The VEV of the dilaton cj) determines the basic string 
coupling constant Qs = e.'^^'^ [!]• One can view this "softening" of all 
coupling constants as a generalization of Einstein's theory of gravity. 
Before Einstein (1915), both the structure of spacetime and the laws 
of local matter interactions were supposed to be "rigid" , i.e. given once 
for all, as absolute structures, independently of the material content of 
the world. Einstein's theory of general relativity introduced the idea 



* Invited talk at the workshop on Varying Fundamental Constants, 3-5 September 
2002 (JEN AM 2002, Porto, Portugal); to appear in a special issue of the Astrophysics 
and Space Science Journal. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



sub.tex; 7/02/2008; 2:11; p.l 



2 



Thibault Damour 



that the structure of spacetime might be "soft", i.e. influenced by its 
material content. However, he postulated ("principle of equivalence") 
that the laws of local physics, and notably the values of all the (dimen- 
sionless) coupling constants {e^ /he, me/mp, . . .), must be kept "rigidly 
fixed" . String theory treats symmetrically spacetime and matter inter- 
actions and suggests that both of them arc "soft". It is amusing to 
note that this generalized correlated "softening" of structures (which 
were traditionally considered as independent and rigid) serendipitously 
shows up even in the mathematical notation used to represent them, 
through a multiple use of the letter "5" : at the tree-level of string theory 
there is a link not only between geometry and gravitation (through the 
unified geometrical field g^,y(a;)), but also between gravitation (g^jy(x)), 
string coupling (gs(a;)), gauge couplings (g(x)) and gravitational cou- 
pling (G(a;)). We refer here to a low-energy Lagrangian density of the 
form {a' = 1) 



Actually such a tree-level Lagrangian (with a massless dilaton (p) is in 
conflict with experimental tests of the equivalence principle. Indeed, 
the dilaton has gravitational-strength couplings to matter which vio- 
late the equivalence principle [2, 3]. For instance, using the results of 
Ref. [3], one derives that the Lagrangian (1) predicts a violation of 
the universality of free fall at the level Aa/a ~ 10~^ (to be compared 
with the present limits 10^^^ [4]), and a time variation of the fine- 
structure constant on cosmological scales: dlne^/dt ~ 10~^°yr~^ 
(to be compared with the Oklo limit ~ 5 x lO^^^yr""^ [see below], or 
with laboratory limits ~ 10~^^yr"^ [5, 6]). 

It is generally assumed that this violent conflict with experimental 
tests of the equivalence principle is avoided because, after supersym- 
metry breaking, the dilaton^ might acquire a (large enough) mass: say 
m^f, > 10~^ eV so that observational deviations from Einstein's gravity 
are quenched on distances larger than a fraction of a millimeter. If that 
were the case, there would also be no possibility to predict any time 
variation of the coupling constants on cosmological scales. There exists, 
however, a mechanism which can naturally reconcile a massless dilaton 
with existing experimental data: this is the deeoupling meehanism of 
Ref. [3] (see also [7]). In the following, we shall review a recent work 
[8, 9] which has extended this mechanism in a manner which comes close 
to reconciling present experimental tests of the equivalence principle 

^ In the following, we use the word "dilaton" to denote the combination of the 
ten-dimensional dilaton and of various moduli which determines the values of the 
four-dimensional coupling constants. 




(1) 
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with the recent results of Webb et al. [10] suggesting that the fine- 
structure constant has varied by ~ 10~^ between redshifts of order 
1 and now. For recent reviews of the flurry of works concerning the 
observational and theoretical aspects of the "variation of constants" 
see [11], and the other contributions to this workshop, notably [12]. 



2. Decoupling mechanism and dilaton runaway 

The basic idea of Rcf . [3] was to exploit the string- loop modifications of 
the (four dimensional) effective low-energy action (we use the signature 
- + ++) 

J \ a a 

-\BF{<t^)F^-...y (2) 

i.e. the 0-dependencc of the various coefficients Bi{(j)), i = g, (p, F, . . . , 
given in the weak-coupling region {e'^ — > 0) by series of the form 

Bi (0) = e~'^ + Cq^ + c^*^ e*^ -|- C2 e^"^ + • • • , coming from the genus 
expansion of string theory. It was shown in [3] that, if there exists 
a special value (pm of (f) which extremizes all the (relevant) coupling 
functions B~^{(fj), the cosmological evolution of the graviton-dilaton- 
matter system naturally drives towards cpm (which is a fixed point 
of the Einstein-dilaton-matter system) . This provides a mechanism for 
fixing a massless dilaton at a value where it decouples from matter 
("Least Coupling Principle"). Refs. [8, 9] considered the case (recently 
suggested in [13]) where the coupling functions, at least in the visi- 
ble sector, have a smooth finite limit for infinite bare string coupling 
gs — 00. In this case, quite generically, we expect 

Bi{ct>) =Ci + 0{e-t>) . (3) 

Under this assumption, the coupling functions are all extremized at 
infinity, i.e. a fixed point of the cosmological evolution is (j)m = -|-oo. [See 
[14] for an exploration of the late-time cosmology of models satisfying 
(3).] It was found that the "decoupling" of such a "run-away" dilaton 
has remarkable features: (i) the residual dilaton couplings at the present 
epoch can be related to the amplitude of density fluctuations generated 
during inflation, and (ii) these residual couplings, while being naturally 
compatible with present experimental data, are predicted to be large 
enough to be detectable by a modest improvement in the precision 
of equivalence principle tests (non universality of the free fall, and. 
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possibly, variation of "constants"). This result contrasts with the case 
of attraction towards a finite value (pm which leads to extremely small 
residual couplings [15]. 

One assumes some primordial inflationary stage driven by the po- 
tential energy of an inflaton field %. Working with the Einstein frame 
metric g^^ = Bg{(f))g^,^, and with the modified dilaton field ip = 
J d(l)[{3/4){B'g/Bg f + B'^/Bg + (1/2) B^/Bg]^/^, one considers an effec- 
tive action of the form 



~ 2 

m 



S = J d^x^ 



2 



(4) 



where nip = l/(47rG) = ACg/a', and where the dilaton dependence of 
the Einstein-frame action is related to its (generic) string-frame depen- 
dence (2) by Fi^p) = B^{^)/Bg{,p) , Vix, ^) = mp^ Bf{<^) y(x, </>). 

Under the basic assumption (3), dip/d4> tends, in the strong-coupling 
limit (p — > +00, to the constant 1/c, with c = (2(7^/(7^)^/^, so that the 
asymptotic behaviour of the bare string coupling is 

gl = e^^e''^. (5) 

Let us consider for simplicity the case where F[ip) = 1 and V{x,'p) = 
X{ip)x^/n with a dilaton-dependent inflaton coupling constant A((^) of 
the form 

X{<p) = X^{l + bxe~'^n, (6) 

where we assume that b\ > 0, i.e that X{(p) reaches a minimum at 
strong-coupling, ip +oo. It is shown in [9] that this simple case 
is representative of rather general cases of (/^-dependent inflationary 
potentials V^x^f)- 

During inflation (ds^ = —dt'^+a^{t) Sij dx^ dx^), it is easily seen that, 
while X slowly rolls down towards x ~ 1, the dilaton ip is monotically 
driven towards large values. The solution of the (classical) slow-roll 
evolution equations leads to 

+ ^ = const. = e^^- + ^ • (7) 

Using the result (7), one can estimate the value i^end of ip at the end of 
inflation by inserting for the initial value Xin of the inflaton the value 
corresponding to the end of self-regenerating inflation [16]. One remarks 
that the latter value can be related to the amplitude dn ^ ^ x 10~^ 
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of density fluctuations, on the scale corresponding to our present hori- 
zon, generated by inflation, through Xin — 5^/n{^H)~^^^^~^'^^ ■ Finally, 
assuming e'^'^^'^ <C e'^'^'^'^'^, one gets the estimate: 

gC'^end ^ i2.5c2 bx {dnr^ . (8) 

A more general study [9] of the run-away of the dilaton during infla- 
tion (including an estimate of the effect of quantum fluctuations) only 
modifies this result by a factor 0(1). It is also found that the present 
value of the dilaton is well approximated by <fend- 



3. Deviations from general relativity induced by a runaway 

dilaton 

Eq. (8) tells us that, within this scenario, the smallness of the present 
matter couplings of the dilaton is quantitatively linked to the smallness 
of the (horizon-scale) cosmological density fluctuations. To be more 
precise, and to study the compatibility with present experimental data, 
one needs to estimate the crucial dimensionless quantity 

aAi^p) = dlnmA{^)/dcp, (9) 

which measures the coupling of (p to a massive particle of type A. The 
definition of aA is such that, at the Newtonian approximation, the inter- 
action potential between particle A and particle B is —Gab ^s/rAB 
where [7, 3] Gab = G{1 + aACXB)- Here, G is the bare gravitational 
coupling constant entering the Einstein- frame action (4), and the term 
aAOLB comes from the additional attractive effect of dilaton exchange. 

Let us first consider the (approximately) composition-independent 
deviations from general relativity, i.e. those that do not essentially 
depend on violations of the equivalence principle. Most composition- 
independent gravitational experiments (in the solar system or in binary 
pulsars) consider the long-range interaction between objects whose 
masses are essentially baryonic (the Sun, planets, neutron stars). As 
argued in [2, 3] the relevant coupling coefficient a a is then approxi- 
mately universal and given by the logarithmic derivative of the QCD 
confinement scale Aqcd(v)) because the mass of hadrons is essentially 
given by a pure number times Aqcd(v')- [We shall consider below the 
small, non-universal, corrections to m^((^) and uaI^) linked to QED 
effects and quark masses.] Remembering from Eq. (2) the fact that, 
in the string frame (where there is a fixed cut-off linked to the string 
mass Ms ^ (a')~^/^) the gauge coupling is dilaton-dependent {g^^ = 
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Bp (if)), we see that (after conformal transformation) the Einstein- 
frame confinement scale has a dilaton-dependence of the form 

Aqcd(v?) ^ B;'/\^) exphSvr^ b^' Bf{v)] Ms , (10) 

where 63 denotes the one-loop (rational) coefficient entering the Renor- 
malization Group running of gp- Here BF{if) denotes the coupling to 
the SU(3) gauge fields. For simplicity, we shall assume that (modulo 
rational coefficients) all gauge fields couple (near the string cut off) to 
the same Bp{ip). Such an assumption is natural in a stringy frame- 
work. Note that we differ here from the line of work of Bekenstein [17], 
recently extended in [18, 19], which assumes that ip couples only to the 
electromagnetic gauge field. The string-inspired assumption of coupling 
to all gauge fields yields the following approximately universal dilaton 
coupling to hadronic matter 

Numerically, the coefficient in front of the R.H.S. of (11) is of order 
40. Consistently with the basic assumption (3), one parametrizes the 
dependence of the gauge coupling gp = Bp^ as 

Bp\p) = Bp\+^)[l-hFe-''^]. (12) 

We finally obtain 

ahad(¥')-406Fce-^'^. (13) 

Inserting the estimate (8) of the value of (p reached because of the 
cosmological evolution, one gets the estimate 

ahad(^end)-3.2^(5^. (14) 

It is plausible to expect that the quantity c (which is a ratio) and 
the ratio hp/hx are both of order unity. This then leads to the nu- 

8 

merical estimate a\^^ ~ 10(5^^^, with 5// ~ 5 x 10~^. An interesting 
aspect of this result is that the expected present value of OiYisA depends 
rather strongly on the value of the exponent n (which entered the 
inflaton potential V{x) x")- In the case n = 2 (i.e. V{x) = \ y^) 
wc have a^^.^^ ~ 2.5 x 10^^, while if n = 4 {y{x) = \ '^X^) we 
have oi\^ ~ 1.8 x 10~^. Both estimates are compatible with present 
(composition-independent) experimental limits on deviations from Ein- 
stein's theory (in the solar system, and in binary pulsars). For instance, 
the "Eddington" parameter 7 — 1 ~ — 2 a^^^ is compatible with the 
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present best limits I7— l|<2xl0 ^ coming from Very Long Baseline 
Intcrferometry measurements of the deflection of radio waves by the 
Sun [20]. 

Let us consider situations where the non-universal couplings of the 
dilaton induce (apparent) violations of the equivalence principle. This 
means considering the composition-dependence of the dilaton coupling 
a A, Eq. (9), i.e. the dependence of a a on the type of matter we consider. 
Two test masses, made respectively of A- and i?-type particles will 
fall in the gravitational field generated by an external mass ruE with 
accelerations differing by 

=2 ■ ~(aA-as)aE- (15) 

\ a J AB aA + aB 

We have seen above that in lowest approximation ua — "had does not 
depend on the composition of A. We need, however, now to retain the 
small composition-dependent effects to aA linked to the (/^-dependence 
of QED and quark contributions to niA- This has been investigated 
in [3] with the result that aA — cthad depends linearly on the baryon 
number B = N + Z, the neutron excess D = N — Z, and the quantity 
E = Z{Z - 1)/{N + Z)V3 linked to nuclear Coulomb effects. Under 
the assumption that the latter dependence is dominant, and using 
the average estimate A{E/M) ~ 2.6 (applicable to mass pairs such 
as (Beryllium, Copper) or (Platinum, Titanium)), one finds that the 
violation of the universality of free fall is approximately given by 

(^) c 5.2 X 10-5 <d - 5.2 X 10-^ (^) ' 5# . (16) 

This result is one of the main predictions of the present model. If one 
inserts the observed density fluctuation 5h ~ 5 x 10-^, one obtains a 
level of violation of the universality of free fall (UFF) due to a run- 
away dilaton which is Aa/a ~ l.3{bF/{b\c))'^ x 10^^^ for n = 2 (i.e. 
for the simplest chaotic inflationary potential V{x) = ^ iti^{4>) x^), and 
Aa/a ~ 0.98{bF/{bxc)f x lO"'' forn = 4 (i.e. for V{x) = iH(l>)x^)- 
The former case is naturally compatible with current tests (at the ~ 
10"^^ level [4]) of the UFF. Values n > 4 of the exponent require 
(within this scenario) that the (unknown) dimensionless combination 
of parameters (bp/ibxc))'^ be signiflcantly smaller than one. 
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4. Cosmological variation of "constants" 

Let us also consider another possible deviation from general relativity 
and the standard model: a possible variation of the coupling constants, 
most notably of the fine structure constant jhc on which the strongest 
limits are available. Consistently with our previous assumptions we 
expect oc Bp^iip) so that, from (12), e^{ip) = e^{^oo) [1 - hper'^'P]. 
The logarithmic variation of (introducing the derivative Lp' = dip /dp 
with respect to the "e-fold" parameter dp = H dt = da/a) is thus given 

by 

c^hpce ^' ~ ^ahadV' • (17) 

The value of (p' depends on the coupling of the dilaton to the two cur- 
rently dominating energy forms in the universe: dark matter (coupling 
ami'f)), and vacuum energy (coupling ay = jd\nV{ip)/d p>). In the 
slow-roll approximation, the cosmological evolution of (p is given by 

(^Ijn + 2^vW = -^mOlm - ^^VOIV , (18) 

where fi^ and J7y are, respectively, the dark-matter- and the vacuum- 
fraction of critical energy density {pc = i3/2)mj,H^). The precise value 
of ip' is model-dependent and can vary (depending under the assump- 
tions one makes) from an exponentially small value {ip' ~ e"'^'^) to a 
value of order unity. In models where either the dilaton is more strongly 
coupled to dark matter than to ordinary matter [21], or /and plays the 
role of quintessence (as suggested in [14]), ip' can be of order unity. 
Assuming just spatial flatness and saturation of the "energy budget" 
by non-rclatistic matter and dilatonic quintessence, one can relate the 
value of ip' = dip/{Hdt) to and to the deceleration parameter 
q = —da/a^: 

ip'^ = I + q - ^n^. (19) 

The supernovae la data [22] give a strict upper bound on the present 
value go- (?o < 0. A generous lower bound on the present value of flm is 
^mo > 0.2 [23]. Inserting these two constraints in Eq.(19) finally yields 
the safe upper bound on the current value of (p' 

(^'o^ < 0.7, i.e. |<^()| < 0.84 . (20) 

On the other hand, Eq. (16) yields the link 

ahad - ±1.4 X 10-S/lOi2 — . (21) 
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Inserting this result in Eq. (17) yields 



dlne^ 
Hdt 



~ ±3.5 X 10 



,-6 




(22) 



which yields, upon integration over p = In a + est = — ln(l + z) + est, 



where = (v(p) ^ 'f^{Po))/{p ^ Po) denotes the average value of if' 

between now and redshift z. If we insert in Eq. (22) the limit Aa/a < 
10~^^, coming from present experimental tests of the universality of free 
fall (UFF) [4], as well as the cosmological constraint (20) on the present 
value of ip', we find that the present variation of the fine-structure 
constant is constrained to be \dlne^/Hdt\ < 3 x 10~^, i.e. \dln.e'^/dt\ < 
2 X 10~^^yr~^. Such a level of variation is comparable to the planned 
sensitivity of currently developed cold-atom clocks [6]. 

However, there are stronger constraints coming from geochemical 
data. Let us first recall that a secure limit on the time variation of 

coming from the Oklo phenomenon is |Ae^/e^| < 10^^ between 
now and ~ 2Gyr ago, i.e. \dlne'^/dt\ < 5 x 10~^^yr~^. This limit 
was obtained in [24] under two very conservative assumptions: 1. A 
conservative interpretation of Oklo data making minimal assumptions 
about the temperature of the reaction zone, and about the possible 
amplitude of variation of the resonance energy of the relevant excited 
state of Samarium 150, and 2. An analysis of the e^-variation of the 
latter resonance energy Er taking into account only the (rather well- 
known) Coulomb effects. We note that Rcf. [25] derived a stronger limit 
on the variation of by replacing assumption 1. above by the non- 
conservative assumption that Er has stayed close to its present value. 
We note also that Ref. [26] derived a stronger limit on the variation of 

by relaxing assumption 2. and by trying to estimate the indirect e^- 
dependence of Er coming from light quark contributions to the nuclear 
binding energy. 

It is important to note that Ref. [26] also derived the strong limit 
lAeVe^l < 3 X 10"^ between now and 4.6 Gyr ago. This limit was 
derived from the Rhenium/Osmium ratio in 4.6 Gyr old iron-rich me- 
teorites by making quite conservative assumptions. In particular, we 
note that, similarly to the Oklo assumption 2. above, this limit uses 
only the Coulomb effects in the /3-decay rate of Rhenium 187. [Ref. [26] 
quotes also stronger limits obtained by trying to estimate the indirect 
e^-dependence of the latter /3-decay rate.] 




(23) 
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Using Eq. (23), the conservative "Oklo" and "Rhenium" Hmits on 
the variation of e^, corresponding to redshifts ^oklo — 0.15 and — 
0.45, can be converted into constraints on the product r {'^')z where r = 
y'lO^^ Aa/a. In round numbers, one finds that these two constraints 
read 



1012 _ |(v7%^o.i5| < 0.2; Wl0i2 — |((^').^o.45| < 0.2. (24) 



As seen on Eq. (18), the cosmological evolution of is driven by two 
quantities: the couphng am of ip to dark matter, and its couphng ay 
to vacuum energy. If one had only the Oklo constraint to cope with, 
one could fine-tune the ratio am/av to satisfy the first constraint (24) 
without constraining the overall magnitudes of a^n and ay ■ This is what 
was done in Ref. [19] (within the different context of Bekenstein-like 
models) to exhibit models satisfying the Oklo and UFF constraints and 
allowing for a variation of around z ~ 1 driven by a large enough Um 
(a la [27]) to explain the observational results of Webb et al. [10]. [Note 
that the implementation of the same idea within the context of dilaton- 
like models [8, 9] leads to a maximal possible variation of which falls 
short, by a factor ~ 4, of the level needed to explain the results of 
[10].] However, the new point we wish to emphasize here is that the 
recently obtained "Rhenium constraint" [26], i.e. the second inequality 
(24), makes it impossible to concoct a fine-tuned ratio am/ctv so as 
to satisfy the two geochemical constraints (24), which correspond to 
quite different redshifts and therefore to significantly different relative 
weights Qrn oc (1 -I- and Cly cx (1 -|- z)°. If we were to consider more 
complicated models, in particular models where the "kinetic term" 
cx (f" must be included in Eq. (18), and allows for an oscillatory 
behaviour (as in local-attractor models [3]), it might be possible to 
fine-tune more parameters (like the phase of oscillation of (p) so as 
to satisfy the two constraints (24). However, such a heavy fine-tuning 
becomes very unnatural. The most natural conclusion is that both am 
and ay must be small enough, so that the quantity ^ylO^'^ Aa/a\{ip')z\ 
is smaller than 0.2 for all redshifts where the cosmologically dominant 
energy forms are dark matter and/or dark energy. Eq. (23) then leads 
to the constraint 



Note that this constraint is about ten times smaller than the claim of 





e^(^) -e^(0)| 



< 0.7 X 10-^ 



ln(l + z) . 



(25) 



[10]. 
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5. Conclusions 

A first conclusion is that the results of Webb et al. [10] cannot be 
naturally explained within any model where the time variation of the 
fine-structure constant is driven by the spacetime variation of a 
very light scalar field tp. On the other hand, the recently explored [8, 9] 
class of dilaton-Iikc models with an attractor "at infinity" in field space 
naturally predicts the existence of small, but not unmeasurably small, 
violations of the equivalence principle. In the case where the dilaton 
if is more significantly coupled to dark matter and/or dark energy 
than to ordinary matter, dilaton-likc models can lead to a cosmological 
variation of as large as Eq. (25). [We recall that this upper bound 
takes into account the two stringent geochemical bounds on Ae^/e^ 
coming from the conservative interpretations of Oklo data [24] and of 
Rhenium decay data [26].] A time variation of the order of the upper 
limit in Eq. (25) might be observable through the comparison of high- 
accuracy cold-atom clocks. It might also be observable in astronomical 
spectral data (if one understands how to explain and subtract the 
systematic effects leading to the current apparent variation of at 
the level AeVe^ -0.7 x 10"^ [10]). 

Finally, an important conclusion of all theoretical models where the 
time variation of is linked to the spacetime variation of a light scalar 
field if (be it the dilaton of string theory [3, 8, 9] or a field constrained 
to couple only to electromagnctism [17, 18, 19]) is that a necessary 
condition for having a fractional variation of larger than 10^® on 
cosmological time scales is to have a violation of the universality of 
free fall (UFF) larger than about 10~^^ (see Eq. (22) in which \ip'\ 
is certainly constrained by cosmological data to be smaller than 1, 
as discussed in [9]). Note that a measurably large violation of the 
UFF is a necessary, but by no means sufficient, condition for having 
a measurably large cosmological variation of e^. Indeed, in Eq. (22) 
or Eq. (23) the value of ip' = d(p/{Hdt) depends on the strength of 
the coupling of to the dominant forms of energy in the universe: 
am and ay - These quantities could well be comparable to the strength 
of the coupling of p to hadronic matter, Ohad; which is contrained 
to be small by UFF experiments (see Eq. (21)). This shows that the 
best experimental probe of an eventual "variation of constants" is to 
probe their spatial variation through high-precision tests of the UFF, 
rather than their (cosmological) time variation (see also [28] for a 
detailed discussion of clock experiments). This gives additional mo- 
tivation for improved tests of the UFF, such as the Centre National 
d'Etudes Spatiales (ONES) mission MICROSCOPE [29] (to fly in 2004; 
planned sensitivity: Aa/a ^ 10~^^), and the National Aeronautics and 
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Space Agency (NASA) and European Space Agency (ESA) mission 
STEP (Satellite Test of the Equivalence Principle; planned sensitivity: 
Aa/a- 10-18) [30]. 
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